Can we hear Kolmogorov spectra? * 
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We study the long-time evolution of waves of a thin elastic plate in the limit of small deformation 
so that modes of oscillations interact weakly. According to the theory of weak turbulence a nonlinear 
wave system evolves in long-time creating a slow transfer of energy from one mode to another. We 
derive this kinetic equation for the spectral transfer in terms of the second order moment. We 
show explicitly that such a non-equilibrium theory describes the approach to an equilibrium wave 
spectrum, and describes also an energy cascade, often called the Kolmogorov-Zakharov spectrum. 
We perform numerical simulations confirming this scenario. 

Introduction.- Since more than forty years it was established that long time statistical properties of a random 
fluctuating wavy system possess a natural asymptotic closure because of the dispersive nature of waves and the 
weakly nonlinear wave interaction 0, Q ■ Indeed this so-called "weak turbulence theory" has shown to be a powerful 
method to study the evolution of nonlinear dispersive wave systems. It results that the longtime dynamics is driven 
by a kinetic equation for the distribution of spectral densities. This method, was developed for surface gravity waves 
0, Q , surface capillary waves Q , plasma waves [f| , nonlinear optics , etc. 

The actual kinetic equation has non-equilibrium properties similar to the usual Boltzmann equation for dilute gases, 
thus it conserves energy, momentum, and exhibits an H-thcorem driving the system to equilibrium, characterized by 
the named Rayleigh- Jeans distribution. Most important, besides the elementary equilibrium (or thermodynamic) 
solution, Zakharov has shown JEj that power law non-equilibrium solutions also arise, namely the Kolmogorov- 
Zakharov (KZ) solutions or KZ spectra, which describe the exchange of conserved quantities (e.g. energy) between 
large and small length scales. 

Experimental evidence of KZ spectra have been found in ocean surface Q and in capillary surfaces waves 0, 0, El • 
On the other h and, numerical simulations of surface waves shown the realization of KZ spectrum for weak turbulent 
capillary waves |ll| and, more recently, for gravity waves |l2j . 

In this article an oscillating thin elastic plate or shell is considered. Adding inertia to the well known (static) theory 
of thin plates one finds ballistic dispersive waves 0] , which interact via nonlinear terms that are weak if the plate 
deformations arc small. A previous work has dealt with solitary wave propagation on the surface of a cilindrical shell 
(l5| due to the balance between dispersive bending waves and nonlinearities. We develop thus a weak turbulence 
theory for the surface deflection. Shortly, the bending waves travel randomly over the system and interact resonantly 
between each other via the weak nonlinearities. The mathematics beyond the resonant condition are formally identical 
to the conservation of energy and momentum in classical gas of particles. In this sense an clastic plate with is formally 
equivalent to a 2D gas of classical particle interacting with a non-trivial scattering cross-section. Indeed, an isolated 
system evolves from a random initial condition to a situation of statistical equilibrium like a gas of classical particles 
does. As usually, in addition to statistical equilibrium for isolated systems, the weak turbulence theory predicts here 
an energy cascade from a source of energy (a driving forcing) to a dissipation scale typically characterized by plastic 
deformations. Moreover, while there is often a lack of direct observations of weak turbulence predictions, we exhibit 
numerically these behaviors for the plate dynamics. 

This dynamics is illustrated in Fig. [J] for an isolated dissipation free system where the plate deformation are shown 
at initial time and after long evolution. 

Theory. - The starting point is the dynamical version of the Foppl-von Karman equations [l^.ll4l| for the amplitude 
of the deformation £ (x, y, t) and for the Airy stress function x( x > V> t)'- 
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FIG. 1: Zoom over a portion of the surface plate deflection ((x,y), the left had image is the initial condition while the right 
hand image represents a late evolution of the elastic plate. The aspect ratio is 1:1:1. 



Where h is the thickness of the elastic sheet. The material has a mass density p, a Young modulus E. Its Poisson ratio 
is a. A = d xx + dyy is the usual Laplacian and the bracket {•, •} is defined by {f,g} = fxx9yy + fyyffxx^^fxygxy, which 
is an exact divergence, so equation preserves the momentum of the center of mass, namely dtt J Ci x i Vi t)dxdy = 0. 
Equation @ for the Airy stress function x(x,y,t) may be seen as the compatibility equation for the in-plane stress 
tensor which follows the dynamics. 

Small plane waves perturbations (£ ~ e i(k-x-u>t) w j^ n x _ (x,y)) of a plane plate are dispersive with the usual 

ballistic behavior of bending waves: ujk = hc\k\ 2 = hck 2 . Where c = J ' i 2 {i-o^ nas dimension of a speed. 

Weak turbulence theory.- Despite the complexity of and J2J) the system presents a hamiltonian structure which 
is straightforward in Fourier space. Defining the Fourier transforms as (k(t) = 57 / ((x,t)e lk x d 2 x (with = C-k)> 
then one gets from J2J): Xk{t) = — 2 ^i {(, (}k where {C,C}fc i s the Fourier transform of {CiC}- The dynamics then 
reads: 
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where d 2 fe2,3,4 = d 2 k 2 d 2 k^d 2 ki and Vi2 ; 34 = 2(2^-)^ ^ fc 1 X |^fi +^fc^ 3 1 !^ fc 4 ^ . The hamiltonian structure becomes evident, if we 
define as canonical variables the deformation £&(£) and the momentum Pk{t) = pdtCk(t). 

The canonical transformation (k = ^|(^U + A*_ k ) an d Pk = — ^ x (Ak — A*_ k ) with X^ = allows to write 

the wave equation in a diagonalized form: + iujkAk = iN^{Ak) where N^{-) is the cubic nonlinear interaction 
term. 

This nonlinear oscillator has two distinct scales of time, the rapid oscillation iujkAk and the weak nonlincarity. 
If the gradient of the deformation is small, the nonlinear terms come up to next order. Then, after the change 
Ak = cik(t)e~ luJkt which removes the rapid linear oscillating term, we obtain: 
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where we define a| with the two possible choices s = +, — corresponding to propagation direction, such that = a k 
while a~ = a*_ k . The interaction term reads: Jk u k*MM = ^ x kiXk 2 X k3 X k4 P 234 Vk 1 ,k2;k 3 M with "P 2 34 the six possible 
permutations between 2, 3 & 4. Next step consists to write a hierarchy of equations for the averaged moments: 
( a fci a fc 2 2 )' { a k a k a k a k^)-> etcetera: (a^a 1 ^) — ie 2 N 3 ((a k \a s k a s k a^)) , etcetera. A multi-scale analysis provides 
a natural asymptotic -over long times- closure, for higher moments: the fast oscillations drive the system close to 
Gaussian statistics and higher moments are written in terms of the second order moment: (afc^^) = rik^ 2 ^ {k\ + k 2 ), 
where rik is called the wave spectrum. 

The wave-spectrum satisfies thus a Boltzmann-type kinetic equation showing the exchange of energy from one mode 
to another in longtime due to the four waves resonance: 
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As for the usual Boltzmann equation, Eq.@ conserves "formally" |2^| the total momentum per unit area P = 
h Jkn k (t)d 2 k and the kinetic energy per unit area £ = hj ujkn k (t) d 2 k and exhibits a if-theorem: let be S(t) = 
J\n(n k )d D k the non-equilibrium entropy, then dS/dt > 0, for increasing time|24j. However, despite the four wave 
interaction type kinetic equation Q, the "wave action" Af = J n k (t)d 2 k is not conserved. The kinetic equation 
PJl describes thus an irreversible evolution of the wave-spectrum towards the Rayleigh-Jeans equilibrium distribution 
which reads, when P = 0: 

^k 

where T is called, by analogy with thermodynamics, the temperature (with units of energy/lenght, that is a force) 
which is naturally related to the initial energy by £q = h J ujkn eq d 2 k = hT J d 2 k. The quantity J d 2 k is the number 
of degrees of freedom per unit surface. Therefore each degree of freedom takes the same energy: hT. Naturally, for 
an infinite system this number diverges (as well as the energy). This classical Rayleigh-Jeans catastrophe is always 
suppressed due to some physical cut-off which corresponds here to dissipation processes for wave-lengths smaller than 
h. Numerical simulations on regular grid provide also a natural cut-off k c = ir/dx with dx the mesh size, which gives 
£q = irhTk 2 for a large plate. 

Kolmogorov Spectra.- In addition, isotropic non-equilibrium distribution solution can also arise They have a 
major importance in the non equilibrium process for the energy sharing among scales. Those solutions can be guessed 
via a dimensional analysis argument but, as it was shown by Zakharov Q, they are exact solutions of the kinetic 
equation. Indeed, integrating over the angles the scattering amplitude | J klk2 k 3 k 4 \ 2 S^ (fei + k 2 + k 3 + £4) a new 
scattering amplitude depends only on the moduli fcj = 
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This integral can be performed providing a resonance condition involving the angle between k 2 & k 3 and ^4 via 
k 2 + k 2 + 2k 2 ■ k 3 = k\ + k 3 + 2k\ki cos ^4. An explicit expression for S is computed and since the degree of 
homogeneity of \ J\ 2 in k is zero, S scales like 1/fc 2 . 

Looking for a power law solution of the form rik = Ak~ a , one has, that the eight terms of the collisional integral in 
the r.h.s. of equation |(2j decomposes into Coll 2 + 2 + Coll 3 +i, where : 
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In Coll 2 + 2 the integration domain is over fl up = {0 < k 2 < k & \Jk 2 — k\ < k 3 < k} and k\ = k\ + k 3 — k 2 while in 
C0II3+1 the integration is over Qdown = {0 < k 2 < k & < k 3 < \Jk 2 — k 2 }, with k\ = k 2 — k\ — k\. 

The collisional terms scaling follow Coll 2+2 = C\(a)k 2 ~ 3a and Colls+i = C 2 (a)k 2 ~ 3a . The coefficients 61/2(0!) are 
pure real functions depending only on a. Both coefficient vanish with double degeneracy at a = 2 indicating that the 
Kolmogorov spectrum coincides with the Rayleigh-Jeans solution JSJ): n^ z ~ i In fact, this degeneracy reveals the 
existence of a logarithmic corrections, similarly to the nonlinear Schrodingcr equation in 2D 0], thus: 
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Here P is the energy flux imposed in the energy cascade between the longwave scales and the short-ones (it has 
dimensions of mass/time 3 ), C and z are pure real numbers. 

For a = and 3a — 4 = the collisional part C0II2+2 also vanishes. Those solutions are the wave action equipartition 
(q = 0) and a second KZ spectrum ~ 1/fc 4 / 3 corresponding to wave action inverse cascade. However, those 
spectrum do not vanish the second part of the collision term C0II3+1, in agreement with the non conservation of 
the wave action mentioned above. Therefore, an important consequence is the non existence here of this second 
inverse cascade ~ 1/fc 4 / 3 , as usually found for four wave interaction systems such as gravity waves or the nonlinear 
Schrddinger equation. Nevertheless, we have observed that for elastic plates the wave action conservation is only 
weakly violated as seen in direct numerical simulations (inset of Fig. ■ 

Numerics.- We are first performing numerical simulations of the full non linear system of PDE Q and J5J with 
no forcing nor dissipation terms. We have taken c = 1, and the ratio h/L, L being the size of a square plate, is 
the only dimcnsionlcss parameter in the numerics. We consider a regular grid with periodic boundary conditions. 
Taking advantage of the structure of the equation in Fourier space a pseudo-spectral scheme is implemented using 
FFT routines namely: (k = —^ICk + {C x}k- The linear part of the dynamics is calculated exactly: £&(£ + At) = 

Cfc (t) cos(wfc At) + sin(cjfe At) . The nonlinear terms in (QJ and are first computed in real space and the integration 
in time is then performed in Fourier space with an Adams-Bashford scheme, which interpolates the nonlinear term of 
d as a polynomial function of time (of order one in the present calculations). Energy is conserved as best as 1/100. 
For initial conditions we have taken: (k ~ Co e ~ k ^ k °e l(pk with ipk a random phase, and a zero velocity field Cfc = 0. 
As time evolves, the random waves oscillates with a disorganized behavior, as shown in FIG. ^ After a long time 
evolution the system build an equilibrium distribution in agreement with the Raylcigh-Jeans nk ~ T/fc 2 spectrum 
which correspond for the plate deflection to: (|Cfc| 2 ) = X%nk = = ph^k 4 - as snown m Fig- 12 
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FIG. 2: Numerical simulation for a square plate with h/L = 10 -3 and using 1024 2 modes. We plot the power spectrum of 
mean deflection (|0c| 2 ) versus wave number k after 1200 time units. The the line plots the Rayleigh- Jeans power law 1/fc 4 . 
The inset plots the magnitude of the violation of wave action conservation. 

Non equilibrium distribution may also be easily observed numerically. One requires to input energy at low wavenum- 
bers (fc < ki n ) and dissipate it at large wavenumbers (fc > k out ) establishing a non- equilibrium Kolmogorov spectrum 
at the window of transparency ki„ < k < k out . This artifact could be easily implemented, adding a Fj. - "fkCk to the 
plate equations. Here the forcing is a nonzero random force for the long waves scales, and 7^ is a fictitious linear 
damping, as in 0. The KZ spectrum is then still observed (see Fig. |3l but it presents a deviation from the 1/fc 4 
spectrum in agreement with a logarithmic correction (inset of Fig. |3J). 

Conclusions.- We have presented a new wave system where the weak turbulence method is successfully applied. 
Numerical simulations exhibit both the convergence towards statistical equilibrium for free systems and the energy 
cascades when forcing and dissipation are introduced. 

As amplitude of the deformations become larger, the elastic plate equations still remain valid, but stretching cannot 
be longer treated as weak perturbation and a "wave breaking" phenomena is expected: energy focuses into localized 
structures as ridges [l9j and conical surfaces (named d- cones) |20j . Amazingly, a regime dominated by ridges shows 
a power spectrum |£fc| 2 ~ 1/fc 4 similar to the weak turbulence spectrum derived here. On the other hand for d- 
cones dominated regimes, as seemingly observed in [2lJ, the expected spectrum should follow |£fc| 2 ~ 1/fc 6 . Finally, 
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FIG. 3: Average power spectrum <|Cfe| 2 ) for the energy cascade, the line plots the power law 1/fc 4 . Inset plots fc 4 (|Cfc| 2 ) vs - 
log k showing a clear logarithmic deviation. 



experimental efforts seem to display a KZ spectrum close to the one predicted by weak turbulence theory [22|. 

In conclusion, since vibrating plates transmit sound to air, we argue that a Kolmogorov spectrum can be heard! 

A. Boudaoud, E. Hamm and L. Mahadevan are acknowledged for communication their experimental results prior 
to publication. This work was supported by FONDECYT 1020359 and 7050143. 



[1] K. Hasselmann, J. Fluid Mech. 12 481 (1962); Ibid. 15 273 (1963). 

[2] D.J. Benney, P.G. Saffman, Proc. Roy. Soc. London A 289, 301 (1966). 

[3] V.E. Zakharov, N.N. Filonenko, Dokl. Akad. Nauk SSSR 170, 1292 (1966). 

[4] V.E. Zakharov, N.N. Filonenko, Zh. Prikl. Mekh. I Tekn. Fiz. 5, 62 (1967) 

[5] V.E. Zakharov, Zh. Eksper. Teoret. Fiz. 51, 686 (1966) 

[6] S. Dyachenko, et al.„ Physica D 57, 96 (1992). 

[7] P.A. Hwang, J. of Phys. Oceanography 30, 2753 (2000). 

[8] W. Wright, R. Budakian, and S. Putterman, Phys. Rev. Lett 76, 4528 (1996). 

[9] E. Henry, P. Alstr0m and M. T. Levinsen, Europhys. Lett., 52, 27 (2000). 
[10] M. Yu. Brazhnikov, G. V. Kolmakov, A. A. Levchenko and L. P. Mezhov-Deglin, Europhys. Lett., 58, 510 (2002). 
[11] A. N. Pushkarev and V. E. Zakharov, Phys. Rev. Lett. 76, 3320 (1996). 

[12] A. I. Dyachenko, A. O. Korotkevich and V. E. Zakharov, Phys. Rev. Lett. 92, 134501(2004). 
[13] L.D. Landau and E.M. Lifshitz, Theory of Elasticity, (Pergamon Press, New York 1959). 

[14] A. Foppl, Vorlesungen iiber technische Mechanik, Bd. 5, Leipzig, 1907, p. 132. and Th. von Karman, Ency. d. math. Wiss., 

Bd. IV. 2, II, Leipzig, 1910, §8. 
[15] J. Wu, J. Wheatley, S. Putterman, and I. Rudnick, Phys. Rev. Lett. 59, 2744 (1987). 
[16] A. Newell, S. Nazarenko, L. Biven, Physica D 152-153 520 (2001). 

[17] V. E. Zakharov, V. S. L'vov and G. Falkovich, Kolmogorov Spectra of Turbulence I (Springer, Berlin, 1992). 
[18] FFTW: www.fftw.org. 

[19] T.A. Witten and H. Li, Europhys. Lett. 23, 51 (1993). 

[20] M. Ben Amar and Y. Pomeau, Proc. Soc. Lond., A 453, 729 (1997); E. Cerda and L. Mahadevan, Phys. Rev. Lett. 80, 

2358 (1998). 
[21] L. Mahadevan, private comunication. 
[22] A. Boudaoud and E. Hamm, private comunication. 

[23] "Formally" means here that the proof requires convergence of any simple integral to the exchange of integration order by 
the Fubinis theorem flr^l . 

[24] For t < one has that dS/dt < 0, driving the system to an equilibrium too, however it should be noticed that kinetic 
equation (|1J is not well posed in the sense that a positive initial condition for rikit — 0) > could become negative. 



